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Abstract 



As shown in arXiv:hep-th/0405282 , the warped deformed conifold has two bosonic massless 
modes, a pseudoscalar and a scalar, that are dual to the phase and the modulus of the baryonic 
condensates in the cascading gauge theory. We reconsider the scalar mode sector, mixing fluctuations 
of the NS-NS 2-form and the metric, and include non-zero 4-d momentum k^. The resulting pair of 
coupled equations produce a discrete spectrum of m\ = —k^ which is interpreted as the spectrum of 
J^'-' = 0^ glueballs in the gauge theory. Similarly, we derive the spectrum of certain pseudoscalar 
glueballs with J^'-'' = , which originate from the decoupled fluctuations of the RR 2-form. We 
argue that each of the massive scalar or pseudoscalar modes we find belongs to a 4-d massive axial 
vector or vector supermultiplet. We also discuss our results in the context of a finite length throat 
embedded into a type IIB flux compactification. 



1 Introduction 



Duality between the cascading SU{k{M + 1)) x SU{kM) gauge theory and type IIB strings on the 
warped deformed conifold [Ij provides a rich yet solvable example of gauge/string correspondence 
O El S] . For earlier work leading up to this duality, see El 13 [5] , and for reviews [10] ■ This 
background demonstrates in a geometrical language such features of the SU{M) supersymmetric 
gluodynamics as color confinement and the breaking of the Z2M chiral R-symmetry down to Z2 via 
gluino condensation [1]. In fact, it has been argued [1] that by reducing the continuous parameter 
QgM one can interpolate between the cascading theory solvable in the supergravity limit and J\f = 1 
supersymmetric SU{M) gauge theory. 

The problem of finding the spectra of bound states at large gsM can be mapped to finding 
normalizable fiuctuations around the supergravity background. This problem is complicated by the 
presence of 3- form and 5-form fiuxes, but some results on the spectra are already available in the 
literature [111 [T2l [T3| [T^ [T5l [TBI I17j . A particularly impressive effort was made by Berg, Haack 
and Miick (BHM) who used a generalized PT ansatz [18j to derive and numerically solve a system 
of seven coupled scalar equations [141 [E] • Each of the resulting glueballs is even under the charge 
conjugation Z2 symmetry preserved by the KS solution (this symmetry was called the X-symmetry 
in [l3]), and therefore has J^'~' = 0"*"+. The present paper will study three other families of glueballs, 
which are odd under the T-symmetry. Two of them originate from a pair of coupled scalar equations, 
generalizing the zero momentum case studied in |TB], and have J^*^ = 0^' . The third, pseudoscalar 
family arises from a decoupled fluctuation of the RR two-form C2 and has J^^ = . 

An important aspect of the low-energy dynamics is that the baryonic U{1)b symmetry is broken 
spontaneously by the condensates of baryonic operators A and B. This phenomenon, anticipated 
in the cascading gauge theory in [HHS], was later demonstrated on the supergravity side where the 
fluctuations corresponding to the pseudoscalar Goldstone boson and its scalar superpartner [13], as 
well as the fermionic superpartner [16], were identified. Furthermore, finite deformations along the 
scalar direction give rise to a continuous family of supergravity solutions [201 [HI 122] dual to the 
baryonic branch, AB = const, of the gauge theory moduli space. 

The main purpose of the present paper is to obtain a deeper understanding of the GHK scalar 
fluctuations [13] and their radial excitations. Our motivation is two-fold. On the one hand, we 
seek an improved understanding of the glueball spectra and their supermultiplet structure. On the 
other, we would like to shed new light on the normal modes of the warped deformed conifold throat 
embedded into a string compactification, which has played a role in models of moduli stabilization 
|23] and D-brane inflation [231 [25] . such inflation models, the reheating of the universe involves 
emission of modes localized near the bottom of the throat, which are dual to glueballs in the gauge 
theory [Ml |271 [25] . 

This paper is structured as follows. In section 2 we construct a generalization of the ansatz 
for the NSNS 2-form and metric perturbations that allows us to study radial excitations of the 
GHK scalar mode. We derive a system of coupled radial equations and determine their spectrum 
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(the details of the numerical treatment are presented in Appendix [D]). In section 3 we show that 
a similar ansatz for the RR 2-form perturbation decouples from the metric giving rise to a single 
decoupled equation for pseudoscalar glueballs. In section 4 we argue that the scalar glueballs we find 
belong to massive axial vector multiplets, and the pseudoscalar glueballs belong to massive vector 
multiplets. Agreement of the corresponding equations is explicitly demonstrated in the large radius 
(KT) limit. In section 5 we give a perturbative treatment of the coupled equations for small mass 
that allows us to study the scalar mass in models where the length of the throat is finite. Review of 
the supergravity equations and of the warped deformed conifold, as well as some technical details, 
are delegated to the Appendices. 

2 Radial Excitations of the GHK scalar 

The ansatz that produced a normalizable scalar mode independent of the four-dimensional coordi- 
nates was [13] 

6B2 = x{r) dg'' , 5Gn = 5^24 = i^{r) . (1) 

Our first goal is to find a generalization of this ansatz that will allow us to study the radial excitations 
of this massless scalar, i.e. the series of modes that exist at non-vanishing A;^ = — m|. Thus, we 
must include the dependence of all fields on x^^. Such an ansatz that decouples from other fields at 
linear order is 

5F^ = ^, 

SB2 = x{x, r) dg^ + d^a{x, r) dx^' A g^ , (2) 
= d6B2 = x'dTA dg^ + d^{x- cr) dx^" A dg^ + d^a' dr A dx^' A g^ , 
5Gi3 = 5G24 = tp{x,T) . 

The ansatz for 6B2 originates from the longitudinal component of a 5-d vector: 

6B2 = {ArdT + A^dxf") A / . (3) 

Requiring the 4-d field strength to vanish, F^^, = 0, restricts to be of the form 0^ acting on a 
function. Then, choosing 

Ar = -X' , A^ = d^{a - x) , (4) 

we recover the ansatz ([2]) up to a gauge transformation. 
Yet another gauge equivalent way of writing ([2|) is 

SB2 = {x-(^) dg^ -a'drAg^. (5) 

The new feature of our ansatz compared to the generalized PT ansatz used in |14t ll5| is the presence 
of the second function in 6B2 which multiplies dr A g^. Terms of this type, which are allowed by 
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the 4-d Lorentz symmetry, turn out to be crucial for studying the modes that are odd under the 
T-symmetry. 

Using S{G-^) = -G-^6GG-^, we find that 6G^'' = SG"^^ = -G^G^^^. The unperturbed 
metric components (see Appendix IB. II for a review of the KS solution) are 

e4/3K(^)sinh2(r/2)/ii/2(r) ' ' 

e4/3K(r)cosh2(r/2)/iV2(^) ' ^ ^ 

^55 _ r^TT _ 6i^(r)^ 

Qf^U ^ ^1/2 ^t.U _ (g) 

In order to find the dynamic equations for the functions x and o" in ([2]) we study the linearized 
supergravity equations below (type IIB SUGRA equations are reviewed in Appendix [A|) . 



2.1 Equations of Motion for NSNS- and RR- Forms 

All the Bianchi identities are automatically satisfied with the ansatz ([2]). Indeed, the relation 
d5H^ = is obvious, and consistent with vanishing dF^ we find that 6H^ A F3 = (using eqs. (|14ip 
and (|145p one can verify that dg^ A F3 = and dr A 5^ A F3 = 0). 
The self-duality equation for F5 reads 

5*^5=0. (10) 

Given that F5 has components along Ag'^ Ag^ Ag'^ Ag^ and along d^xAdr, our adopted deformation 
of the metric does not affect to first order. 

Even though the variations of the forms -F3 and are zero, the deformations of their Hodge 
duals (5 * -F3 and 6*F^ will in general be non-zero because of the deformations of metric components. 
In the equation for F^ 

d5 * F3 = F5 A dHs , (11) 

the product F^ A SH^ vanishes identically. From the explicit form of F^ we see that the Hodge dual 
of the first two terms in (jl4ip will be a closed form 6 * F^ = A{x, r) d'^x A dr A (...). The third 
term in (jl4ip . F'dr A {g^ A g^ + g'^ A g^), is not affected by the deformation of the metric, and thus 
d6 * F3 = is satisfied identically. 

The remaining equations are nontrivial. In particular 

d6*H3 = 0, (12) 

turns out to be more complicated than the equation for F^. The variation 

6*H3 = *dH3 + Sg * H3 (13) 
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consists of two parts: *6H3 accounting for the deformation of the form itself, and 5g*H^ arising 
from the deformation of the Hodge star. ExpUcit calculation shows that 



5g * 



{G^^f IG^^I df,a' *4 dxf" A dg^ A dg^ , 



-G G^^ [f'G^^ + d^x A dg^ A g^ . 



(14) 

The four-dimensional Hodge star *4 is taken w.r.t. the standard Minkowski metric. Differentiating 
this expression for 6 * and equating to zero the coefficients multiplying linearly independent forms 
gives three equations: 



d'^x A dg^ A dg^ 



11 /~«33 



2CG 



G^^ /i2 UW , 



(15) 



d^x A dr A dg^ A g^ : dr 



-GG^^G^'^ G^^ 



gsMa' 



[f'G'' + k'G'']^P + x' 



+ 



*4(ix^ A dr A dg^ A dg^ 



G^^ G^^ Oiix - a) = , (16) 
2^/^G^^G^^h'2^^{x-a) + ^r{V^{G^^fh^^^a'j=0, (17) 



where we have substituted for the warp factor IG^'^j = /12 (no summation over is implied). Not 
all of these equations are independent. Indeed, using ([15]) equation ([16]) simplifies to 



V {\/^ (G^^)2 /l3 □4CJ'} + 2 /il □4(;^ - fj) = . 



(18) 



d. 



This is exactly what we obtain by acting on (|17p with 5^ and contracting indices. Thus only (|15p 
and ()17p are independent. The coefficient functions in these equations are given by (we have dropped 
some inessential constant factor in \/—G): 

4K(r)2 



2 (sinh 2r - 2r) 

16 

e8/3/i(r)K(r)2sinh2T ' 
6if(r)^ 



^11 ^33 



£4/3 

4 



Taking into account these expressions, equations ([15]) and ([17]) read 

2{gsMa' 



-7^%^ + X' = ^ e'/' h{r) K{rf sinh^ r ^4^' , 

£4/3 y^/l(T) 16 



9m(x - c^) + ^ K{Tf dr[K^ sinh^ r a^a'} = . 



(19) 
(20) 

(21) 
(22) 
(23) 

(24) 
(25) 
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2.2 Einstein Equations 

The first order perturbation of the Ricci curvature tensor is given by 



SRij — — (^ — 6Ga°'.ij — 6Gij-^a'^ + 6Gai;j'^ + SGajii"") , 



(26) 



where covariant derivatives and contractions of indices are performed using the unperturbed metric. 
The first term in this expression vanishes because the metric perturbation is traceless. The remaining 
three terms combine to give the only non-zero perturbations 5Ri3 = 5R2a'- 

1 



5^13 = -^i^ sinh(T)z 



K" Ih" z" (K'f 

1 1 V — 

K 2 h z 



K' z' fh' K' z' , 

2—- + cothr - + 4— + 2- ) + 2 
K z \ h K z 



K' h' 

2 /i2 +T77^ 

1 4 1 



sinh(T)2 9 sinh(r)2/i:6 



hi{T)Ksmh{T) UiZ 



p4/3 



K sinh Tz 



(Ksinh(r))^ (In/i)' 
(Ksinh(r))2 



+ 



{Ks\nh{T)Y z' 
[K sinh(r))2z 



^ 4 cosh(T) 



sinh(r)2 <d sin^rf 3iv:3sinh(r) 



—h{T)K sinh r \I\4z . 



(27) 



where z{x,t) is defined by 

V'(x,r) = h^/^Ksmh{T)z{x,T) = 2-^/^[smli{2T) - 2T]^/^h^/^ z{x,t) . 



(28) 



The source terms on the right hand side of the Einstein equation Rij = Tij ()120p are due to 
the deformations of the metric and B2 form. It turns out that the only nontrivial deformations are 
those with indices 13 or 24, with dTi^ = ST24. Say, for the 13 component STi^ we have the following 
contributions: 



1 
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\G^^ Hi2t 5H^2t + G^^ SHi4^r -f^34r] G' 

]{gsMa')G'' [G^' f + G^^ k'] x' 



SGiFiabcdPs"''') = ^ [G^^f {G'^f G^' (F12345)' , 
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SciHiabHs"''') 



2 [Hi35 SG^^ G^^ + Hi2r H^iT ^G^'^ G^^] 
1 

2' 



; [(i/135)' - H^2r i^34r] G^^ G^^ ^ 



(A; - fY - f'k' 



gii ^33 ^55 ^ ^ 



(29) 
(30) 



(31) 
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f 5G(i^la6i^3'^') = f [i^l25 i^345 .^G^^ + (^G^^ G^^] 



5^ 



[(Fi3.)' - i^l25 i^345] G^^G^^ G^^V 

{gsMa'f - F(l - F)] G" G^^ G^^ ^ ^ 



(32) 



+i G^^ {k - ff + (G^i)2 + (G33)2 (1 - Ff 



Denoting 



5T^ 



13 



_2(^11 (^33 i7'2]^. 



^i(r)+^2(r) V'(a;,r) +S(r)x'(x,r) 



where Ai stands for the contribution from F5, we get 

3(5,Ma')^ (rcothr - l)2[sinh(2T) - 2r]4/3 



B{t) 



sinh6(T) 



2l/3g20/3/i5/2 

3 {gsMa' f 
8e4/i3/2 sinh^ r L' 

-3 (a.Ma- '^-'^"'■W 



6^/3 /i(t) sinh^ r 
Then ehminating x' with the help of psp yields 

3 (gsMa')'' (r coth r - 1)^ [sinh 2r - 2t]^/3 



5Ti 



13 



22/3 g20/3/j2 
X 



sinh^ r 



(33) 
(34) 

(35) 

(36) 
(37) 



3 (gsMa'f (sinh(2T) - 2t)1/3 

z{t) + — T7T7 ^ ' ^ a. . X 



3 cosh4r - 8t sinh 2r - 8r^ cosh 2t - 8 cosh 2t + IGr^ + 5 



• 21/3 e^h 



sinh'^(r) 



cosh(4r) + 8(1 + r^) cosh(2T) - 24t sinh(2T) + 16r^ - 9 

□40-'(x,r) 



9 QsMa' sinh2r - 2r 



16 e4/3 



sinhr 



£4/3 



JC'^ sinh r 



1 (/i')^ 1 h" K' h' ^ h' 



9 gsMa' sinh2T - 2r 
16 £4/3 sinhr 



(38) 
□4a' . 



As mentioned above the perturbations (5Ti3 = 5T2i are the only non-zero components of 5Tij. 
Equating (j27p and (j38p we obtain the final form of the linearized Einstein equation. 
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2.3 Two Coupled Scalars 

Combining the equations for the field strengths and the Einstein equations we have the system 

{gsMa') '"'^]t^^.~J^ ^ + X' = ^ e'/' Kr) K{rf sinh^ r U,a' , (39) 
e^l-^ smh^ r 16 

- ^) = -\K{Tf dr[K^ sinh^ r d^a'] , (40) 

((i^sinhr)^z') g4/3/j / 2 8 1 4 coshr \ 

(Ksinhr)2 6i^2 VsinhV OKSsinhV 3/^3sinhVy 

3 , sinh 2r — 2t „ , , 

+- {gsMa') —2 U^a' . (41) 

16 smh T 

Note that x can be ehminated between ([39]) and (fiOll . Further, a change of variables 



5 = zE:sinh(r) , (42) 

4/3 

-u; = _ sinh(r) V , (43) 



gsMa 

leads to a more symmetric pair of equations 

,, 2 e^/^/i^ 3(c/^Ma')^sinh2r-2r^ 

cosher + 1 e^/^/i^ 8sinh2r-2r , 

IL" ^ T-;rlAAW = 5 Z. 45 

sinhV 6^2 9 K^sinh^T ^ ' 

Introducing the dimensionless mass-squared m2 according to 

we can rewrite the equations for z and w as 

„ cosh^ r + 1 9 /(r) 16 , , 

This is a system of coupled equations which defines the mass spectrum of certain scalar glueballs 
with positive 4-d parity. The natural charge conjugation symmetry of the KS background is the 
T-symmetry, under which these modes are odd. Therefore, we assign J^^ = 0"* to this family of 
glueballs 

In the massless case these equations lead to the GHK solution [13j. If we assume 04 = —^;^ = 
m\ = 0, then there are two solutions |l3j, zi = cothr and 22 = Tcothr — 1. The solution for z 



^ For comparison, the glueballs found in |141ll5j are 0'*'^. The glueballs whose spectrum comes from the minimal scalar 
equation [11] resulting from the analysis of graviton fluctuations are 2++. The axial vector U{l)ii fluctuations [17] give 
rise to 1^^ glueballs whose masses are also determined by the minimal scalar equation. 
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which is non-singular at the origin is z = rcothr — 1. Substituting it into the second equation, we 
find 

coshV + 1, 16 , , , 22/3 8 ,^ ^ . , 

w" 2 w = — K{t) (t coth r - 1) = /'(r) sinh r . (49) 

sinh r 9 9 

The two solutions of the homogeneous equation are wi = 1/ sinhr and W2 = (sinh2T — 2T)/sinhT; 

both of them are singular either at zero or at infinity. This means that the regular solution of the 

inhomogeneous equation is uniquely fixed. With the Wronskian W {wi,W2) = wiw'2 — w'iW2 = 4, we 

can find a general solution 

w{t) = - ^{«)i(r) [Ci - J^dx r{x) sinhx] + ^r) [C2 + f dx I'{x) sinhx] } . 

(50) 

Integrating by parts and choosing the particular homogeneous solution to make w well behaved at 
both zero and infinity we get 

, . 2^/3 8 1 



W[T] 



9 sinh T 



dx I {x) sinh^ X . (51) 



Alternatively, if we start with equation (3.18) in the GHK paper [13], which determines /2, and 
introduce w = /2(r)Er2 sinhr, (jl9|) . up to a rescaling of the right hand side. Thus, we can 

use the solution (3.25) of [13] for /2(t) to read off the result ([5T|) . 



Let us also note that the non-zero w in the zero momentum case A:^ = is not in contradiction 
with the GHK solution. This is because w enters ([2|) only through d^a which is zero as long as the 
momentum vanishes. 



2.4 Numerical Analysis 

To determine the spectrum of glueballs in the field theory, we need to solve the eigenvalue problem 
for in the infinite throat limit. This system of equations (j47p . ()48p does not seem amenable 
to analytical solution and we employ a numerical approach to find the spectrum of normalizable 
solutions. It is convenient to use the determinant method, which generalizes the standard shooting 
technique to a system of several coupled equations (see [H]). The detailed description of the 
numerical analysis as well as of the subtleties specific to the system ()47p . (j48p is given in Appendix iDl 
The result is that the spectrum consists of two distinct series, each with a quadratic growth of 
for large n. These series are interpreted as the radial excitation spectra of two different particles. 
The lowest eigenvalues {fh? < 100) for these spectra are shown in Table [TJ The quadratic fit for 
spectrum I is 

mf„ = 2.31 + 1.91 n-F 0.294 (52) 
For spectrum II (we drop the lowest eigenvalue when doing the fit) 

mfj„ = 0.36 + 0.14 n + 0.279 n^. (53) 
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Table 1: Non-zero eigenvalues with fh? < 100. There 
by quadratic polynomials in the eigenvalue number n i 

Spectrum I 



n 


ml 


n 


ml 


n 


ml 


n 


ml 


1 


4.53 


5 


19.1 


9 


43.3 


13 


76.9 


2 


7.30 


6 


24.4 


10 


50.8 


14 


86.7 


3 


10.7 


7 


30.1 


11 


58.9 


15 


97.1 


4 


14.6 


8 


36.4 


12 


67.6 







m~ 




5 10 15 



e the two distinct spectra. Both spectra can be fitted 
e red line in the plots). 

Spectrum II 



n 


ml 


n 


ml 


n 


ml 


n 


ml 


1 


0.129 


6 


8.06 


11 


30.1 


16 


65.1 


2 


0.703 


7 


11.2 


12 


35.5 


17 


73.9 


3 


1.76 


8 


15.0 


13 


42.1 


18 


83.3 


4 


3.33 


9 


19.3 


14 


49.2 


19 


93.3 


5 


5.43 


10 


24.1 


15 


56.9 
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It is interesting to compare these results with those found for the 0^^ modes by Berg, Haack and 
Miick (BHM) [15]. The conventions of [15j correspond to a particular choice of the KS parameters 
(see Appendix IB. 3p . and the relation between the masses is 

m%HM = (3/2)2/3/(0) 0.9409 . (54) 

Using this relation one can convert the mass eigenvalues to the BHM normalization. We note that 
the lightest glueball we find, the first entry from spectrum II in Table 1, has w-^/^jvf ~ 0.121. 
For comparison, the lightest O"*"^ eigenvalue found in [15] has itl^bhm ^ 0.185. The fact that 
the 0"' sector has the lightest glueballs may be qualitatively understood as follows. Roughly 
speaking, glueball masses increase with the dimensions of the operators that create them. The 
lowest dimension operator from the 0"''+ sector is the gluino bilinear TrAA of dimension 3, but the 
0^ sector contains an operator of dimension 2, namely Tr(^^ — BB). 
Converting the asymptotics of the two spectra to BHM units, we find 

"1/ BHM ~ 2.17 + 1.79 n + 0.277 (55) 
"1?/ Bi?M~ 0.34 + 0.13 n + 0.262 (56) 
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The coefficients of the quadratic terms are close to those found in [15]. The quadratic dependence 
on n, which is characteristic of Kaluza-Klein theory, is a special feature of strongly coupled gauge 
theories that have weakly curved gravity duals (see [29] for a discussion). Note that m| is obtained 
from through multiplying by a factor ~ Ts/{gsM), where Tg is the confining string tension. 
Thus, for n <^ ^fg^M these modes are much lighter than the string tension scale, and therefore 
much lighter than all glueballs with spin > 2. Such anomalously light bound states appear to be 
special to gauge theories that stay very strongly coupled in the UV, such as the cascading gauge 
theory; they do not appear in asymptotically free gauge theories. Therefore, the anomalously light 
glueballs could perhaps be used as a 'special signature' of gauge theories with gravity duals if they 
are realized in nature. 

One may be puzzled why the spectrum in Tabled] does not include the GHK massless mode. This 
is because in solving the coupled equations (I47p . ()48p we required that both wave- functions z and w 
vanish as r ^ cxd. This excludes the GHK zero mode which grows as i ~ r. On the other hand, this 
growth is a lot slower than the exponential growth found for generic solutions. The meaning of the 
GHK mode as the baryonic branch modulus seems to be well established since even the solutions 
at finite distance along this modulus are available [201 121j . Thus, the GHK scalar zero-mode should 
be normalizable with a proper definition of norm. In fact, the GHK pseudoscalar and its fermionic 
superpartner are normalizable \yi\ I16j: therefore, the supersymmetry of the problem implies that 
the GHK scalar is normalizable as well and is part of the spectrum. 



3 Pseudoscalar Modes from the RR Sector 

The type of ansatz used in section O works even more simply for the RR 2-form field: 

5F5=0, 

= x{x-, t) dg^ + d^a{x, r) dx^ A g^ , 
(5F3 = d6C2 = x'dT A dg^ + d^{x - cr) dxf" A dg^ + d^a' dr A dx^" A g^ . 



(57) 



This ansatz is similar to, but somewhat simpler than the GHK pseudoscalar ansatz [13] which 
involved mixing with 5F^. Since dF^ A = 0, now it is consistent to set 5F^ = 0. We also have 
-F5 A SFj, = 0, so it is consistent to take SH^ = 0. Finally, one needs to study mixing with metric 
fluctuations. At a first glance it seems that 5Gi2 and (5G34 might need to be turned on, but a more 
detailed analysis shows that their sources vanish: 

5ri2=Fi3.5F23- + 5Fi4.F24^ = ^G^^G^s^fV-FV] = 0, (58) 
<5T34 = F3i,(5F4i" + 5F32.F42- = 0. (59) 
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Thus, the perturbation (j57p decouples from all other modes, and the only non-trivial linearized 
equation is 

d*5F3 = 0. (60) 
The calculation we need to perform is the same as in section 2.1, except we now set 1^ = and find 

= ^ ^(^) ^(^)' si^h' ^ °4Ct' , (61) 
16 

df,{x-(T) + ^K{Tfdr[KUmh^Td^a']=0. (62) 
Eliminating x and changing variables as before, 

4/3 

* = ^rr7^'sinMT)V, (63) 

we find 

cosher + 1 e'^/^/j 

-r^ w + -TTT^O^w = . 64 

smh^ T DA ^ 

Again, after introducing the dimensionless mass as in (j46p . we get a non-minimal scalar equation 

„ cosh^ T + 1 9 /(r) , , 

w" 5 w + m^-^-J^w = 0. (65) 

sinh T A(rj"= 

Since the 4-d parity operation includes sign reversal of RR fields, we identify the family of glueballs 
coming from this eigenvalue problem as pseudoscalars whose J^*^ quantum numbers are . 

If we set m = the solution regular at small r is (sinh 2r — 2r) / sinh r. Since this blows up at 
large r we conclude that this equation does not contain a massless glueball. A simple numerical 
analysis using the shooting method allows one to find the mass spectrum. The lowest eigenvalues 
(m^ < 100) are listed in Table [2j The quadratic fit is 

mf^j^ = 0.996 + 1.15 n + 0.289 ; (66) 

in the BHM normalization it is given by 

mjii BHM = 0-938 + 1.08 n + 0.272 . (67) 

The spectrum can be reproduced with good accuracy using a semiclassical (WKB) approxi- 
mation. The effective potential in (j65p is singular at r = which does not allow us to use the 
conventional WKB approximation. Yet we can cast the equation (j65p in the form Q1Q2W = m'^w, 
where Qi are first-order differential operators and then consider an equation Q2Q1W = m?w, which 
must give rise to the same spectrum up to a zero mode. Namely, in our case this means that for A 
such that 

,9 cosh^ T + 1 , , 

A^ + A' = ——^ , (68) 

smh T 
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Table 2: Non-zero eigenvalues with < 100 in the RR sector. This spectrum can also be fitted by a quadratic 
polynomial (red line). 



80 - 



n 


ml 


n 


ml 


n 


ml 


n 


ml 


1 


2.41 


5 


14.0 


9 


34.8 


13 


64.7 


2 


4.47 


6 


18.0 


10 


41.4 


14 


73.7 


3 


7.08 


7 


23.2 


11 


48.6 


15 


83.2 


4 


10.3 


8 


28.7 


12 


56.4 


16 


93.3 



equation (j65p shares the spectrum with an equation 



I(t) 

w'' - (B'^ + B')w + m'^—^w = 0, (69) 



A general solution of (j68p reads 



2sinh^r , 

^ = -cothT + — -. 71 

cosh T smh r — r + C 

For (j69p to be non-singular at the origin C has to be non-zero. For a finite C the potential is regular 
everywhere but not monotonic and (j69p admits a zero mode. A most convenient choice is to take 
infinite C, which reduces A to A = — cothr. In this case the WKB approximation is applicable 
in it simplest form (see [11] for similar considerations) and yields the same result as the shooting 
method up to the third digit. 



4 Organizing the Modes into Supermultiplets 

The pseudoscalar Goldstone mode and the massless scalar found in [13j belong to a 4-dimensional 
chiral multiplet. These fields appear as the phase and the modulus of the baryonic order parameters 
that vary along the baryonic branch. When a long KS throat is embedded into a Calabi-Yau 
compactification with fluxes, the baryonic U (1) symmetry becomes gauged and a super symmetric 
version of the Higgs mechanism is expected to take place. The axial vector U{1)b gauge field 'eats' 
the pseudoscalar mode and acquires a mass degenerate with the mass of a scalar Higgs. These fields 
constitute the bosonic content of a massive M = 1 axial vector supermultiplet. 
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In the present paper we explicitly constructed the massive modes that are radial excitations of 
the GHK scalar. It is, of course, interesting to find the supermultiplets they belong to. We will argue 
that each of these scalar radial excitations is also a member of a massive axial vector supermultiplet. 
Similarly, each pseudoscalar glueball found in section [3] is a member of a massive vector multiplet. 
To prove these facts we would need to demonstrate the existence of the J^*-" = 1"' glueballs 
degenerate with the 0"' glueballs found in section [2l as well as of 1 glueballs degenerate with 
the glueballs found in section [3l Unfortunately, constructing decoupled equations for vector 
supergravity fluctuations around the KS background is a difficult task. Instead, we will provide 
some evidence for our claims by studying axial vector and vector fluctuation equations in the large 
radius (KT) limit (setting a' = g<j = 1, = and M = 2; see Appendix IB. 2p . 

First we reconsider the simple decoupled pseudoscalar equation from the RR sector ([64|) and 
argue that its superpartner is given by the four-dimensional vector Ai in 

5B2 = Al^g\ (72) 

where we have chosen the ansatz so that the corresponding radial component ^ dr A vanishes. 
The equation for d * H-^ implies (with primes denoting derivatives with respect to r) 



,3 



TT *4 A 





3 



4:T hv 

-TT — Ci4 *4 ^4 ^1 = , (73) 

6 D 



and *4 Ai = 0, i.e. the vector is divergence-free. Since the Laplacian acting on such a vector is 
□4 = — *4 (^4 *4 (^4 (note the Minkowski signature of the four dimensional metric), we find 

A[ Ai + 04^1 = . (74) 



6 



3 6 



Defining a new variable Ai = rAi, it is easy to see that its equation of motion. 



- a; -Ai + '-^UiAi = Q , (75) 



3 



3 



coincides with the KT-limit of the equation for the decoupled pseudoscalar w, once we identify 
r ~ e^l'^e^l'^ . In fact, if we make the same ansatz (172p in the full KS background, the equation of 
motion for A\ = sinh tAi resulting from the terms ~ d^x Adr ALO2ALO2 in d* = is precisely 
as in ([M]) : 

d^ ~ cosher + 1 ~ 6^/3/1^ ~ ^ 

However, this ansatz is not closed in the KS case. The Bianchi identity for is not satisfied, so 
this NSNS vector must mix with RR excitations of F3 and/or in the KS background. It would 
be interesting to solve this mixing problem. 

Let us now turn to the massive axial vector superpartners of the coupled scalars (f44l) , (f45|) found 
above. We make the following ansatz, which is similar to the one studied in |30j . 

SC4 = Bi AL03 + F2 AuJ2 + Ki Adr Auj2, (77) 
5C2 = CiAg^ + D2 + EiAdr, (78) 
6B2 =H2 + JiAdr; (79) 
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where Bi,Ci, Ei, Ji, Ki are axial vectors and D2, F2, H2 are two- forms in four dimensions. We 
choose to spht the six degrees of freedom residing in the two-form into a vector and a dual vector, 
e.g. D2 = d4{. . .) + *4d4^Di. The degrees of freedom contained in the former (exact) part are in fact 
the same as those in Ei, so we can simply write D2 = *40?4-Di without loss of generality. Similarly, 
F2 = *4d4Fi and H2 = *4,diHi, and the corresponding exact parts can be absorbed into the vectors 
Ki and Ji, respectively. 

The equations of motion then imply that Bi and Ci have to be divergence-free: ^4 *4 Bi = 
(i4 *4 Ci = 0. If this were not the case their divergences would simply couple to additional scalars 
6C4, ~ dr A iL>3 and 6C2 drAg^, respectively, but we will not consider this here (i.e. as for Ai above 
we choose as gauge in which these radial components vanish) . In fact we will assume that all vectors 
in our ansatz are divergence-free, and that the terms appearing in the RR- and NSNS-potential^ 
are eigenstates of the Laplacian ^4 = — *4 d^ *4 ^4 with eigenvalue m^. 

We have relegated the details of the derivation of the equations of motion to Appendix [Cl 
Splitting the equations obtained from (jll9p into exact and coexact parts w.r.t. the four-dimensional 



derivative operator d^ shows that the vectors Ei , Hi and Ki decoupL 
the remaining vectors read 



hr 



3 ^ 



/ 3 

+ - D^Bi 
r 

hr ^ ^ 



□4D1 



Ji + -Cu 

r 



The resulting equations for 

(80) 



6 ^ 
54 1 



3 D 



-□4F1 



54 



-F{ --Bi 

r 



3 log — Ji 



hr^ 
'5A 



+ 31og— 



Fi + -Bi 
r 



□4F1 



hr 



B[ = —UiJi + Z\og-U^Di 
54 r* 



(81) 
(82) 
(83) 
(84) 
(85) 



where (|8T|) . (|83|) and (j84j) hold modulo terms annihilated by d^. It is easy to see that (j84j) and (|85|) 
imply ([5D|) . so the latter is not independent. We thus have the five coupled equations for the five 
vectors Bi,Ci, Di, Fi and Ji . 



^I.e. we demand that for example □4-B1 = m^Bi, but for the vectors derived from two-forms, such as we only 
impose the weaker condition ^4 *4 d^Di = m? *4 diDi. 

^More precisely, we set Ei = —3 log ^ Hi = r log ^ Ki , and find a single second order differential equation obeyed 
by these fields. Thus we have found another decoupled vector, but this is not the one we are looking for. Given this relation 
between them, Ei, Hi and Ki do not mix with the other vectors. 
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In the massless case, our ansatz includes the pseudoscalar found in [13j. Putting 



Ci = -f2{r) dia{x) 
OiDi = fi dia{x) , 

F, = J,=0 



(86) 
(87) 



/i/irlog — d4a{x) , 



(89) 



for some constant /i and a four-dimensional massless pseudoscalar a(x), all equations of motion are 
satisfied provided 

' T"^ "1 ' ^T" Q T" 

(90) 



4r 

y/2 



9 _ r 
--/ilog — 



in perfect agreement with the literature. 

Now we would like to consider massive excitations however, and find axial vector-like solutions 
to the equations (I80p - (j85p which give rise to the superpartners of the massive scalar excitations of 
(jM]) and (P5|) . In particular, changing variables to Wi = rCi equation ([82]) becomes 



hr"^ 2 6 , 

Wi + — □4VF1 = -04^1 -T-^B[. 
9 r n,r^ 



(91) 



Thus we can identify Wi with w in ()45p , which suggests setting the right hand side of this equation 
proportional to the counterpart of z/r. Hence we define 



Zi = UiFi -j-B[. 

nr 

Using (|80|) and (i8T|) one can deduce that this new field obeys 



-Z[\ + — 04^1 = - 04^^! + - 04 (Ji + D[) 



(92) 



(93) 



Our reduced ansatz containing five axial vectors is still too general. In order to match the 
spectrum of the scalar particles found above, we need to impose an additional constraint to reduce 
the number of dynamical vectors obeying independent second order differential equations to two. 
The correct constraint for our purposes is given by 



□4(Ji + I)'i) =-U4Wi. 



(94) 



In order to show that we can consistently impose this relation we need to examine the remaining 
equations. First of all, with this constraint ()85p reads 



^1 



hr^ 



hr^ 



-— U^D^ + — 04^^! + Slog -U4D1 
54 18 



Adding ([83j) and ([8^ . using the constraint and the fact that Wi is a mass eigenstate we find 



-Ts" 



+ ^\og—Wi + 



2Jo 



h r 
~5A 



□4^1 



0. 



(95) 



(96) 
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Eliminating □4-D1 between the last two equations we obtain a second order differential equation 
containing only Wi and Zi. A non-trivial fact is that this equation is identical to (j9ip . This 
relies heavily on the precise expression for the warp factor (jlSOp . and shows the consistency of the 
constraint equation with the equations of motion. 

Finally, introducing a symbol for the other combination of the vectors Fi and Bi that appears 
in the equations of motion 



Yi = Fi + -Bi, 



equation (j8^ implies 



□4^1 



z[ --n^Di. 

r 



In summary, we have the two coupled dynamical equations 



r 
3 



3 ^ 



hr^ 2 
— 04^1 = - □4VF1 
9 r 



r 
3 



Wi + — UiWi = -Zi , 



(97) 
(98) 

(99) 
(100) 



which determine Wi and Zi. In terms of these □41)1 is determined by ()96]) . Ji by (j94l) . and □4!! 
by (f98|) . Equations ([99]) . (jlOOp are precisely the KT hmit of the scalar equations ([M]), (|i5]) up to a 
rescaling of the fields by a numerical factor]^ 

Since the KT limits of their equations of motion agree, we thus argue that the axial vectors 
Zi and Wi are the superpartners of the coupled scalars z and w found above, and their massive 
excitations combine into vector multiplets. 



5 Effects of Compactification 

Now we will embed the KS throat into a flux compactification, along the lines of |31], and estimate 
the mass of the Higgs scalar. Generally, glueballs are dual to the normalizable modes localized near 
the bottom of the throat, and one does not expect them to be strongly affected by the bulk of the 
Calabi-Yau. This is indeed the case for all the massive radial excitations found in sections [2] and 
[3j We will see, however, that the case of the GHK scalar is more subtle and exhibits some UV 
sensitivity. 

To model a compactification, WG will introduce b UV cut-off on tlie radial coordinate, Tmax* 

We 

also need to include a deformation of the KS solution introduced by bulk effects. On the field theory 
side this corresponds to perturbing the Lagrangian of the cascading gauge theory by some irrelevant 

^ Looking at (|44p one might have expected the term 2z/ sinh^ t to give rise to a term proportional to Zi/r^ in ((99)) . 
but in fact this is not the case because it is too small to be seen in the KT hmit. In the KS background it arises from a 
subleading term in the variation of the Ricci tensor (P7)) but such terms that are asymptotically suppressed by powers of 
r compared to the leading terms are not taken into account in the KT metric. Indeed, if we write ansatz ^ in the KT 
background and follow the same strategy as we did for the full KS background, the term proportional to z/r^ does not 
appear in the Einstein equations. 
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operators. Here we are not interested in classifying all of them but rather model the compactification 
effects in the simplest way by considering one perturbation which simulates the main features of the 
compactified solution. We consider a shift of the warp factor 5h = const which corresponds to the 
dimension 8 operator on the field theory side |32[ [55] . This also has a simple geometrical meaning: 
the warp factor of the compactified solution is a finite constant in the bulk of the Calabi-Yau and 
therefore should not drop below a certain value along the throat. 

Let us introduce a small parameter 6 which shifts the rescaled warp factor, /(r) I{t) + 5, and 
consider the system ()47p - ()48p in perturbation theory near fh? = 0: 

z = zq + m^zi , (101) 
w = Wo + fh^wi , (102) 
Zq = rcothr — 1 , (103) 

22/3 g ]^ f' 

wo(t) = — / dx I (x) s'mh^ X . (104) 

9 smh r Jq 

At leading order in m? 

zi = (r cothr — 1) / dxu(a;) cothx — cothr / dx u{x) {x coih. x — 1) , (105) 
Jo Jo 

1 r, / ,sinh2x-2x sinh2T-2r f°° , , , I , , 

wi = —- — / dxv[x) / dxvyx) — - — . (106) 

4 sinhr /n sinhx 4 sinhr L. sinhx 





I{t) 9 5 

<^'^=- -]o{t) + 4:2^ ^^^^ " ' 

v{r) = -0 + f K{r) h-^,wo. (108) 

Keeping in mind that for large t, u —2' one finds the asymptotic behavior 

zi(r)~ -2-2/35 Tdx (r - x)a;e2-/3 ^ ^^^^2r/3_ ^^gg) 

Jo 4:2^/-^ 

This yields v ~ —2'^/^6Te'^/^ and 

1 r , , ,sinh2x-22; sinh2r-2r /■°° , , , 1 9 2^/35 

wi = --——— / dxvo{x) — — — / dxvo[x) . ~ — - — re^/''.(110) 

4 smhr Jo smhx 4 smhr J^. smh 2; 8 

Finally, up the first order in the mass squared and 5: 



This suggests that for generic boundary conditions the cut-off value 

Tmax^ -log5=^/2m^ (112) 

This prediction can be tested numerically. In order to do this one can specify some small rh and 
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log m 




10 15 20 25 30 35 



Figure 1: The dependence of logm on Tmax is linear with the slope equal to -1/3. The three lines shown 
correspond to S = I, 6 = 0.01 and S = 0.0001. 

plot the determinant 

detf'^^") ''^^A , (113) 
Yu;i(r) W2{t) J 

of the two linearly independent solutions regular at r = as a function of r. The first zero marks 
the point Tmax such that there is a regular solution with z(Tmax) = '"^(''"max) = 0. Hence Tmax is the 
corresponding cut-off value. As Fig. [1] shows, the relation ()112p holds for Tmax large enough that 

is small. 

Let us consider a simple model of compactification where the throat is embedded into an asymp- 
totically conical space that terminates at some large cut-off value Tmax- To calculate the mass from 
(|114p we need to know S as well as Tmax- The former is the asymptotic value of the (rescaled) warp 
factor. The point where the field theory warp factor approaches S marks the UV cutoff of the field 
theory 

I{tuv) ~ rc/ye-'^^^^/^ ~ S . (115) 

Using this in (jll4p we find m? ~ e'^^'^'^^~^'^™'"=)/^. This shows that the Higgs mass becomes para- 
metrically small only for Tmax 3> 2tjjv This is not satisfied in general; the geometry requires only 
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that Tmax > because tjjv is the length of the throat embedded into a CY space. With the ratio 
between the UV and IR scales of the field theory around 4 • 10^ |24j we estimate that Tjjy ~ 25 
[2T]. The cut-off 

''"max can be related to the warped volume of the Calabi-Yau which, in a singular 
conifold approximation, is 

C = Vol(ri'i) drh{r)J-^^ , (116) 

Jo V det 5^1,1 

where r ~ e^l'^e^l'^ . The integral from zero to ruy is the warped volume of the throat, and from 
ryy to Tmax is the bulk volume. Assuming that the latter dominates, 

- ^^e'l\9sMo!f [rLx " r%y\ . (117) 

Requiring Tmax ^ 50 leads to an enormous ^ far larger than, for example, ~ 5^a''^ in [24j . 

Thus, while for Tmax ^ 2r[/y the Higgs scalar becomes parametrically lighter than the other 
normal modes, in compactifications with realistic parameters it may actually be heavier. This is 
due to the special feature of its wave function z which grows linearly with r in the throat. The only 
conclusion we can draw from our simplified model of compactification is that this mode is rather 
UV sensitive, so to determine its mass we need to know the details of the compactification. 
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A The Type IIB Supergravity Equations 

Here we succinctly list the equations of motion required to study RR and NSNS 2-form perturbations. 
Since the dilaton and RR scalar do not enter at linear order, we set them to zero. 
Bianchi identities: 

di^3 = , 

ARz = ^, (118) 

dFg = H3AF3. 
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Dynamic equations: 

d*H3 = -g^^F5AF3, 

d*F, = F,AHs, (119) 

Einstein equation: 

„2 1 1 2 2 

Tj rp "s TP T7I abed I ZT TJ oh r< u Tjabc \ "s rp tj ah "s rp rpabc /-i on^ 

ilij — J-ij — ■^i'iabcd^j + -^J^iabHj — -^^ijJ^abcJ^ + -^^iab^j — -^'-^ij-t'abc-t' ■ 

B Review of Warped Deformed Conifolds 
B.l The KS Solution 

The ten dimensional metric for the KS solution is 

dsfo = h{T)-^/^{-dt^ + dx"^ + dy^ + dz^) + /i(r)^/2(is^ , (121) 

where 



dSa 



3K^ 

(122) 

is the usual warped deformed conifold metric. The volume form is 

vol = g^h^^^ sinh^ rdt A dx^ A dx"^ A dx^ A dr A A g"^ A A A . (123) 

The one-forms are given in terms of angular coordinates as 

— o — 



1 ^ ^ 2 
5 = ^ , 9 



V2 ' ^ x/2 ' 



o + 4 + 

5 = 7^ , 9 - 



where 



V2 ' " ^/2 ' 

/ = (124) 



= — sin^id^i , = d6i , 

= cosipsm.62d(j)2 — s\n.ipd02 , 
= sin -0 sin ^2(^<?!'2 + cos V"^^2 > 
6^ = ^-0 + COS 9id(f>i + cos 6'2#2 • (125) 
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Also 



and 



Note that 



The warp factor is 



where 



dg' 



, 1 COSt/l 1 9 

sin^i 

de^ = , 

, o 4 c COS 6*1 4 1 

sm ^1 

,4 COS^l o 1 

de^ = -e^ A -i A e\ 

smt'i 

de^ = -e^ A + A e"^ 



1 



V2 



cot 01 (5^+5^) a/, 



(5' - 5^) a/ --^ cot 01 5^ a/, 
~{g'-g')Ag'-^ cot 0i (5^ + 5') A / , 

-(5^a/ + 5^A5'). 



(sinh(2r) - 2r)V=^ 



2V3sinhr 

Mr) = ig.Mar2'/'e-'/'l{r) , 



°° , X coth X — 1 , . , , ^ X 1 /■I 

o?x 2 (sinh(2x) - 2a;) . 

sinh X 



The NSNS two-form field and corresponding field strength are 

B2 = '-^mr)g'Ag' + k{r)g'Ag'], 



(126) 
(127) 
(128) 

(129) 
(130) 

(131) 
(132) 

(133) 

(134) 
(135) 

(136) 
(137) 
(138) 

(139) 



H. = dBo 



9sMa' 



dr A ifg' A g^ + k'g^ a g^) + -{k - f )g^ A {g' A g' + g' A g^) 



while the RR three-form field strength is 

Ma' 



F3 = 



2 

Ma' 



g'Ag'Ag^ + d[F{T){g' A 5' + 5' A /)] 



g"" A g"" A g\l - F) + g'' A g' A g^F + F'dT A {g' A g^ + g^ A g'] 



, (140) 



(141) 
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The auxiliary functions in these forms are 

^ , , sinh T — T 
Ft) 



Some useful identities are 



2 sinh r 

/(r) = "7^:"\ coshr-l), (142) 
2 smh T 

, , , r coth r — 1 

^l''") = — TT^i (cosh r + 1) . 

2 smh r 

A;-/ = 2F', 

/' = (1-F) tanh2(r/2), (143) 
/c' = F coth2(T/2). 

The five-form field strength is given by 

Fs = (1 + *) ^2 A F3 . (144) 

We also note that 

dg' = -{9'/\9' + 9'/\9') , (145) 

and 

dg^ A dc/^ = -2g^ A g"^ A A . (146) 



B.2 The KT Solution 

The KT solution [8] corresponds to the large r limit of the more general KS solution. For simplicity 
we take gs = a' = 1, M = 2 and = 0. In terms of the radial coordinate r ~ g^/Sgr/s ^j^g j^'p 
background is given by 



ds^ = (-dt^ + dx"^) + \/h{dr'^ + r^ds^n) , 



iJs = - dr A W2 
r 



i?2 = 3 log — UJ2 



F5 = (l + *)52 AF3 = Slog- 



r* L 



CJ2 A UJz 



-^3 — '^S ) 

54 



d X Adr 



(147) 
(148) 
(149) 



The warp factor is given by 



h{r) 



81 

8^ 



1 + 4 log ■ 



(150) 



and the conifold metric is 



dsl.. = hg'f + W^g 



i\2 



i=l 



(151) 
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The volume form is given by 



, Vhr^ ,4 
vol = — - — a X f\ L02 f\ u)^ f\ dr . 
54 

Here we have introduced the two harmonic forms, 

^2 = ^(5^ Ag^ +3^ a/) = ^(sin6'i(i6li A #1 - sin 6l2d6'2 A (i^2) , 



(152) 

(153) 
(154) 



B.3 BHM Normalization 

Here we show how to find the conversion factor between the dimensionless mass squared fh? and 
the mass in the normalization of Berg, Haack and Miick [15j. We note that the BHM conventions 
correspond to the KS solution with an extra relation between e and M. The authors of [15] use the 
notations of the general PT ansatz (as given in Eq. (3.8) of [14]): 



ds^ = e^P-'^dsl + (e^+f + a^e^-<^){el + e^) + e^-s[ei + - 2a(eie3 + 6264)] + e" 



6p-x 2 
65 ' 



After settin 



a = tanhy 



1 



cosh T 



, e ^ = cosh y = coth r ; 



(155) 
(156) 

(157) 



it reduces to the KS form 



ds^ =e^^+^P-^r]ijdx'dx^ + 



sinhr 



coth r(ei+ 62 + 61 + 64) + 



sinhr 



(6163 + 6264) + e' 



-&p—2x 



{dr^ + el) 



The radial KS coordinate r is introduced according to 



= e-^Pdr 



Note that in the KS notation the conifold metric (jl22p can be rewritten as 



dsl 



^^(dr^ + (65)^) + ^ cosh r(6i + 63 + 6^ + 64) + 6163 + 6264 



(158) 
(159) 

(160) 



In terms of r, the PT variables necessary to describe the metric for the KS background solution 
take the form 



e^ = tanh(T/2), 

-e*^P+2^ = cothT ^ 
3 



sinh^ r ' 



^2xy3-4p^g~2/3^2g<I.ojj 



r sin 



h^/^r. 



(161) 
(162) 
(163) 
(164) 



^The Papadopoulos-Tseytlin [TB] variables are {x,p,y,^,b,hi,h2)- 
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In the BHM normalization 

.2/3 I{t 



These equations give for the coefficients 

2 

Comparing these coefficients with those of the KS solution we finc|§ 

A/3 

i - 0-1/3 1>o/2 J 

6*0/2 _ i 

"2 ■ 

This yields e'^/^/M'^ = Then using (|46|) we get for the four-dimensional mass in 

BHM normalization 

6 e^/^ 

^%HM = "^4 = ^ ^ = (3/2)2/^ lo ■ 

C Equations of Motion for Vector Superpartners 

With the ansatz (j77p . the deformations of the field strengths are 
5H^ = - *4 U/^Hi + [*/^d4^H[ + ^4 Ji) A dr , 

*(5F3 = — — □4-f/"i A LJ2 A A dr + — - {d4^H[ - ^id^Ji) A A W3 ; 

54 54 ' 

6F3 = d^Ci A - C[ Adr A - Ci A dg^ + {*4dAD[ + d^Ei) A dr + d^ *a d^Di , 

*(5F3 = *4 d4Ci A c<j2 A W2 A dr H *4 A 0^2 A a;2 + - *4 Ci A dg^ Ag^ Adr 

6 6 3 

+ — — (ci4^'i - *4C^4-E'i) A u;2 A a;3 — — □4-D1 A u;2 A 0^3 A dr ; 

54 ^ ' 54 



5F. = + *6J^' 



5 1 



^JTg = (d4Bi - B[a dr) AL03+ {- *4 □4F1 + *4d4Fi' A dr) AL02+ d^Ki Adr A 

3 2> hr 

*6J^5 = - *4 d4Bi A UJ2 A dr + — *4 B[ Auj2 O4F1 Auj^ Adr 

r hr 3 

+^ {d^Fi - *4d4Ki) ACU3. 



UJ2 



^We set gs = a' = I according to [15] . 
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The equations of motion that result from this ansatz are as follows. The Bianchi identity for 
gives 

r 



*4 □4-/^! = — *4 □4-H'l 



*4 □4-61 



hr 



*4 B[ 



*4 □4D1 



hr 



d^n^Fi + 



■ {d4^F[ - *4d4,Ki) = *4d4,H[ + diJi + - d^Ci . 



From the equations for *F^ in (jll9p we have 



hv^ T 

*4 □4£'i = —3 log — *4 □4i7i 

54 



—— *4 □4C1 


— ((i4l)l-*4d4^l) 



' 4r 



9 



— *4 C*! + -:r *4 C*! = *4 □4-?^l + T-^ *4 -S^ , 

+ — - d4n4-Di = *4(i4/i:i - d4Fi' 



54 



-3 log — (*4(i4i/i + (i4 Ji) d4i?i , 



and from the equations for */73 



(179) 
(180) 
(181) 

(182) 
(183) 

(184) 



d4n4-f^i + 



— *4 0471 = -3 log — *4 □41)1 + *4n4Fi - — *4 ^l' , (185) 

54 r* nr 

fxT^ -1 / J. 

— (<i4i?( - *4d4 Ji) = 3 log — {*4d4D[ + (i4^i) 

- 041: J 

- *4 diF[ - diKi - - *4 d^Bi . (186) 
r 



D Numerical Analysis: Finding the Spectra 

A standard method of finding the spectrum of a single second-order differential equation is the 
shooting technique. For a system of several coupled linear equations the shooting method has to be 
generalized [15]. Here we will focus on the subtleties specific to the system of equations (H?!) and 
l8]) . The idea of the calculation (called the determinant method [15]) is to set the initial conditions 



ziir) \ / Z2(r) 

at infinity corresponding to the two solutions regular at infinity, | j and I j , and 

wi{t) ) \W2{'r) 

extend them numerically to small r. Then the matrix 

%(0) %(0)\ (j^^, 

becomes degenerate at the critical points (eigenvalues) in the spectral parameter space. 
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Let us find the asymptotic behavior of regular and singular solutions near both zero and infinity. 
At small r equations ([17|) and ([18|) decouple, 

z"-4^ = 0, (188) 
w"-^w = 0. (189) 

There are the two regular solutions with z, w ^ t'^ and the two singular solutions with z, w ^ l/r. 
For large r we have 

z" = rh'j^e'^/'w, (190) 
w"-w= ^ e-^'^~z. (191) 



The asymptotic behavior of the two regular solutions is 






(z^ 








l = ( 




\W2) 





and the singular solutions are 



-2^/^ (r-f)e- 



r/3 





(~zA 




) 1 













64-21/3 « 

e"^ 



16-21/3 « 



(192) 



(193) 



A particular subtlety of this setup is that at large r the two singular solutions don't diverge 
equally fast: one of them grows exponentially while the other is only linear in r. This makes it 
difficult to start shooting from zero: imposing the regularity condition at infinity would require 
vanishing of both linear and exponential terms. To cancel the linear term in the presence of the 
exponential one is difficult to do numerically. That is why for this particular system it is convenient 
to start shooting from large r, since both singular solutions at zero share the same behavior (~ l/r). 
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